Structure of A = 10 — 13 nuclei with two- plus three-nucleon interactions from chiral 

effective field theory 
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Properties of finite nuclei are evaluated with two-nucleon (NN) and three-nucleon (NNN) in- 
teractions derived within chiral effective field theory (EFT). The nuclear Hamiltonian is fixed by 
properties of the A — 2 system, except for two low-energy constants (LECs) that parameterize the 
short range NNN interaction. We constrain those two LECs by a fit to the A = 3 system bind- 
ing energy and investigate sensitivity of 4 He, 6 Li, 10 ' 11 B and 12 ' 13 C properties to the variation of 
the constrained LECs. We identify a preferred choice that gives globally the best description. We 
demonstrate that the NNN interaction terms significantly improve the binding energies and spectra 
of mid-p-shell nuclei not just with the preferred choice of the LECs but even within a wide range 
of the constrained LECs. At the same time, we find that a very high quality description of these 
nuclei requires further improvements to the chiral Hamiltonian. 

PACS numbers: 21.60.Cs, 21.30.-x, 21.30.Fe, 27.20.+n 



The nuclear strong interaction has proven to be compli- 
cated and replete with ambiguities. However, chiral per- 
turbation theory (ChPT) Q provides a promising bridge 
to the underlying theory, QCD, that could remove ambi- 
guities. Beginning with the pionic or the nucleon-pion 
system Q one works consistently with systems of in- 
creasing nucleon number @, 0, One makes use of 
spontaneous breaking of chiral symmetry to systemati- 
cally expand the strong interaction in terms of a generic 
small momentum and takes the explicit breaking of chiral 
symmetry into account by expanding in the pion mass. 
Thereby, the NN interaction, the NNN interaction and 
also 7rN scattering are related to each other. At the 
same time, the pion mass dependence of the interaction 
is known, which will enable a connection to lattice QCD 
calculations in the future Nuclear interactions arc 
non-perturbative, because diagrams with purely nuclc- 
onic intermediate states are enhanced Therefore, the 
ChPT expansion is performed for the potential. Solving 
the Schrodinger equation for this potential then automat- 
ically sums diagrams with purely nucleonic intermediate 
states to all orders. Up to the fourth- or next-to-next-to- 
ncxt-to-leading order (N 3 LO) of the ChPT, all the LECs 
can be determined by the A = 2 data with the exception 
of two LECs that must be fitted to properties of A > 2 
systems. The resulting Hamiltonian predicts all other 
nuclear properties, including those of heavier nuclei. We 
demonstrate that this reductive program works to pre- 
dict the properties on mid-p— shell nuclei with increasing 
accuracy when the NNN interaction is included. 

We adopt the potentials of ChPT at the orders 
presently available, the NN at N 3 LO of Ref. and the 
NNN interaction at N 2 LO @]. Since the NN interaction 
is non-local, the ab initio no-core shell model (NCSM) 
0, U| is the only approach currently available to 
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FIG. 1: Relations between Co and ce for which the binding 
energy of 3 H (8.482 MeV) and 3 He (7.718 MeV) are repro- 
duced, (a) 4 He ground-state energy along the averaged curve, 
(b) 4 He charge radius r c along the averaged curve. Dotted 
lines represent the r c uncertainty due to the uncertainties in 
the proton charge radius. 



solve the resulting many-body Schrodinger equation for 
mid-p-shell nuclei. In this paper, we use the NCSM to 
evaluate binding energies, spectra and other observables 
for 6 Li, 10 ' n B and 12,13 C. We also present our results for 
the s-shell nuclei 3 H, 3 He and 4 He. We use the A = 3 
binding energies to constrain the two unknown LECs of 
the NNN contact terms, cr> and Ce [12] ■ We then inves- 
tigate sensitivity of the A > 3 nuclei properties to the 
variation of the constrained LECs. Our approach differs 
in two aspects from the first NCSM application of the 
chiral NN+NNN interactions in Ref. [12| which presents 
a detailed investigation of Li. First, we introduce a reg- 
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FIG. 2: Dependence on the cd with the ce constrained by 
the A = 3 binding energy fit for different basis sizes for: 



°Li quadrupole moment, 1U B B(E2;3+0 



l+0)/B(E2;3^0 
1+0) ratio, and the 12 C B(Ml;0+0"-> 1+1). The HO fre- 
quency of Kl — 13, 14, 15 MeV was employed for 6 Li, 10 B, 
12 C, respectively. In the inset of the 12 C figure, the conver- 
gence of the B(Ml;0 + — * 1 + 1) is presented for calculations 
with (using cd = —1) and without the NNN interaction. 



ulator depending on the momentum transfer in the NNN 
terms which results in a local chiral NNN interaction. 
Second, we do not use exclusively the He binding en- 
ergy as the second constraint on the cd and ce LECs. 

The NCSM casts the diagonalization of the infinite di- 
mensional many-body Hamiltonian matrix as a finite ma- 
trix problem in a harmonic oscillator (HO) basis with an 
equivalent "effective Hamiltonian" derived from the orig- 
inal Hamiltonian. The finite matrix problem is defined by 
-ZVmax, the maximum number of oscillator quanta shared 
by all nucleons above the lowest configuration. We solve 
for the effective Hamiltonian by approximating it as a 
3-body interaction [H EO] based on our chosen chiral 
NN+NNN interaction. With this "cluster approxima- 
tion", convergence is guaranteed with increasing -/V max . 

It is important to note that our NCSM results through 
A = 4 arc fully converged in that they are independent 
of the cutoff, N max , and the HO energy, ft£l. For heavier 
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FIG. 3: 1X B excitation spectra as function of the basis space 
size iV max at hQ = 15 MeV and comparison with experiment. 
The isospin of the states depicted is T=l/2. 



systems, we characterize the approach to convergence by 
the dependence of results on iV max and Ml. 

Fig. Q] shows the trajectories of the two LECs cd — ce 
that are determined from fitting the binding energies of 
the A = 3 systems. Separate curves are shown for 3 H and 
3 He fits, as well as their average. There are two points 
where the binding of 4 He is reproduced exactly. We ob- 
serve, however, that in the whole investigated range of 
cd — ce, the calculated 4 He binding energy is within a 
few hundred keV of experiment. Consequently, the de- 
termination of the LECs in this way is likely not very 
stringent. We therefore investigate the sensitivity of the 
p-shell nuclear properties to the choice of the cd — ce 
LECs. First, we maintain the A = 3 binding energy con- 
straint. Second, we limit ourselves to the cd values in 
the vicinity of the point ~ 1 since the values close to 
the point cd ~ 10 overestimate the 4 He radius. 

While most of the p-shell nuclear properties, e.g. exci- 
tation spectra, are not very sensitive to variations of cjj 
in the range from -3 to +2 that we explored, we identi- 
fied several observables that do demonstrate strong de- 
pendence on ce- In Fig. [2] we display the 6 Li quadrupole 
moment that changes sign depending on the choice of cd, 
the ratio of the B(E2) transitions from the 10 B ground 
state to the first and the second l + state, and the 12 C 
B(M1) transition from the ground state to the 1 + 1 state. 
The B(M1) transition inset illustrates the importance 
of the NNN interaction in reproducing the experimen- 
tal value 0. The 10 B B(E2) ratio in particular changes 
by several orders of magnitude depending on the cu vari- 
ation. This is due to the fact that the structure of the 
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two 1+0 states is exchanged depending onco. Using ex- 
trapolation, we can see that the best overall description 
is obtained around the Co « — 1. This observation is also 
supported by excitation energy calculations as well as 
by calculations of other transitions. We therefore select 
cd = —1 and, from Fig. [TJ ce = —0.346 for our further 
investigation. 



We present in Fig. [3] the excitation spectra of n B as 
a function of N niax for both the chiral NN+NNN, (top 
panel) as well as with the chiral NN interaction alone 
(bottom panel). In both cases, the convergence with in- 
creasing 7V max is quite good especially for the lowest-lying 
states. Similar convergence rates arc obtained for our 
other p— shell nuclei. 
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FIG. 4: States dominated by p-shell configurations for 10 B, 
n B, 12 C, and 13 C calculated at jV max = 6 using Ml = 15 MeV 
(14 MeV for 10 B). Most of the eigenstates are isospin T—0 or 
1/2, the isospin label is explicitly shown only for states with 
T—l or 3/2. The excitation energy scales are in MeV. 



We display in Fig. [4]the natural parity excitation spec- 
tra of four nuclei in the middle of the p— shell with both 
the NN and the NN+NNN effective interactions from 
ChPT. The results shown arc obtained in the largest 
basis spaces achieved to date for these nuclei with the 
NNN interactions, 7V raax = 6 (6M7). Overall, the NNN 
interaction contributes significantly to improve theory 
in comparison with experiment. This is especially well- 
demonstrated in the odd mass nuclei for the lowest few 
excited states. The celebrated case of the ground state 
spin of 10 B and its sensitivity to the presence of the NNN 
interaction is clearly evident. There is an initial indica- 
tion in these spectra that the chiral NNN interaction is 
"over-correcting" the inadequacies of the NN interaction 
since, e.g. l + and the 4+0 states in 12 C are not only in- 
terchanged but they are also spread apart more than the 
experimentally observed separation. While these results 
display a favorable trend with the addition of NNN in- 
teraction, there is room for additional improvement and 
we discuss the possibilities below. 

These results required substantial computer resources. 
A typical iV max = 6 spectrum shown in Fig. 3] and a 



set of additional experimental observables, takes 4 hours 
on 3500 processors of the LLNL's Thunder machine. We 
present only an illustrative subset of our results here. 

Tabic Q] contains selected experimental and theoretical 
results for 6 Li and A = 10 — 13. A total of 71 experi- 
mental data are summarized in this tabic including the 
excitation energies of 28 states encapsulated in the rms 
energy deviations. Note that the only case of an increase 
in the rms energy deviation with inclusion of NNN inter- 
action is 13 C and it arises due to the upward shift of the 
| state seen in Fig. |4l an indication of an overly strong 
correction arising from the chiral NNN interaction. How- 
ever, the experimental | may have significant intruder 
components and is not well-matched with our state. 

We demonstrated here that the chiral NNN interaction 
makes substantial contributions to improving the spectra 
and other observables. However, there is room for further 
improvement in comparison with experiment. We stress 
that we used a strength of the 27r-exchange piece of the 
NNN interaction, which is consistent with the NN inter- 
action that we employed. Since this strength is some- 
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TABLE I: Selected properties of 6 Li, 10 ' n B and 12 ' 13 C from 
experiment and theory. E2 transitions are in e 2 fm 4 and Ml 
transitions are in The rms deviations of excited state 

energies are quoted for the states shown in Fig.|4]whose spin- 
parity assignments are well established and that are known 
to be dominated by p-shell configurations. The total energy 
rms is for the 28 excited states from Fig. [4] Results were 
obtained in the basis spaces with A max = 6 (8 for 6 Li) and 
HO frequency hVl = 15 MeV (13 MeV for 6 Li J 14 MeV for 
10 BV The experimental values are from Ref. [jj, HE EE EE 

EE EE IE Hi- 



Nucleus /property 


Expt. 


NN+NNN 


NN 


e Li : |£(1+0)| [MeV] 
Q(l+0) [e fm"] 


31.995 


32.63 


28.98 


-0.082(2) 


-0.124 


-0.052 


fi(l 1 0) [fi N \ 


n coo 


+U.OOD 


+U.840 


E (1 + Cl) fMeVl 


2.186 


2.471 


2.874 


B(E2;3+0 1+0) 
B(E2;2+0 1+0) 
B(M1;0+1 -> 1+0) 


10.69(84) 
4.40(2.27) 
15.43(32) 


3.685 
3.847 
15.038 


4.512 
4.624 
15.089 


B(M1;2+1 -> 1+0) 


0.149(27) 


0.075 


0.031 


iu B : \E(3+0)\ [MeV] 
r v [fm] 


64.751 
2.30(12) 


64.78 
2.197 


56.11 
2.256 


Q(3+0) [e fm 2 ] 


+8.472(56) 


+6.327 


+6.803 


At(3+0) [fx N ] 
rms(Exp - Th) [MeV] 
B(E2;l+0 3+0) 
B(E2;l+0 -> 3+0) 


+1.801 

4.13(6) 
1.71(0.26) 


+1.837 
0.823 
3.047 
0.504 


+1.853 
1.482 
4.380 
0.082 


B(GT;3+0 -» 2+1) 
B(GT;3+0 2+1) 


0.083(3) 
0.95(13) 


0.070 
1.222 


0.102 
1.487 


U B: |S(|7i)| [MeV] 


76.205 


77.52 


67.29 


rpdi |) [fm] 


2.24(12) 


2.127 


2.196 


Q(|-i) [efm 2 ] 


+4.065(26) 


+3.065 


+2.989 


rms(Exp - Th) [MeV] 


+2.689 
- 


+2.063 
1.067 


+2.597 
1.765 


B ( E2 'li 5 li D 


2.6(4) 


1.476 


0.750 


B(GT;2 1 2 21 2) 


0.345(8) 


0.235 


0.663 


B ( GT ;2i 2^212) 


0.440(22) 


0.461 


0.841 


B(GT,2 1 2 ' 21 s) 


0.526(27) 


0.526 


0.394 


B(GT; 5l 2 22 2) 


0.525(27) 


0.762 


0.574 


T2((~ ir p. 3 1 5 1 \ 

V ' '21 2 " 29 2 ' 


0.461(23) 


0.829 


0.236 


i2 C : |£(0+0)| [MeV] 
r p [fm] 


92.162 
2.35(2) 


95.57 
2.172 


84.76 
2.229 


Q(2+0) [e fm 2 ] 
rms(Exp - Th) [MeV] 
B(E2;2+0 -> 0+0) 


+6(3) 
7.59(42) 


+4.318 
1.058 
4.252 


+4.931 
1.318 
5.483 


B(Ml;l+0 -> 0+0) 


0.0145(21) 


0.006 


0.003 


B(M1;1+1 -> 0+0) 


0.951(20) 


0.913 


0.353 


B(E2;2+1 -» 0+0) 


0.65(13) 


0.451 


0.301 


13 C: |S(|-i)| [MeV] 


97.108 


103.23 


90.31 


r P {~ 21 |) M 


2.29(3) 


2.135 


2.195 


rms(Exp - Th) [MeV] 


+0.702 


+0.394 
2.144 


+0.862 
2.089 


B ( E2 ;fi 2 ^ 5i 1) 


6.4(15) 


2.659 


4.584 


B(M1;| L | - | L |) 


0.70(7) 


0.702 


1.148 


B (GT;2 1 2^212) 


0.20(2) 


0.095 


0.328 


B (GT;2 1 2 ~* 21 2) 


1.06(8) 


1.503 


2.155 


B(GT;2 1 2 -* 22 2) 


0.16(1) 


0.733 


0.263 


B(GT; 2 - 1 2 ^ 22 2) 


0.39(3) 


1.050 


0.221 




0.19(2) 


0.400 


0.151 


Total energy rms [MeV] 




1.314 


1.671 



what uncertain (see e.g. Ref. [l2|), it will be important 
to study the sensitivity of our results with respect to this 
strength. Further on, it will be interesting to incorpo- 
rate sub-leading NNN interactions and also four-nucleon 
interactions, which are also order N 3 LO (22)]. Finally, we 
plan to extend the basis spaces to N max = 8 (8SO) for 
A > 6 to further improve convergence. 

Our overall conclusion is that these results provide ma- 
jor impetus for the full program of deriving the nuclcon- 
nucleon interaction and its multi-nuclcon partners in the 
consistent approach provided by chiral effective field the- 
ory. It is straightforward, but challenging, to extend this 
research thrust in the directions indicated. The favor- 
able results to date and the need for addressing funda- 
mental symmetries of strongly interacting systems with 
enhanced predictive power, firmly motivates this path. 
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